Non exponential spin relaxation in magnetic field in quantum wells with random 

spin-orbit coupling 



M.M. Glazoyi and E.Ya. Sherman^ 

'A. F. loffe Physico- Technical Institute, Russian Academy of Sciences, 194021 St. Petersburg, Russia 
'^Department of Physics and Institute for Optical Sciences, University of Toronto, 
60 St. George Street, Toronto, Ontario, Canada M5S 1A7 

We investigate the spin dynamics of electrons in quantum wells where the Rashba type of spin-orbit 
coupling is present in the form of random nanosize domains. We study the effect of magnetic field 
on the spin relaxation in these systems and show that the spatial randomness of spin-orbit coupling 
limits the minimum relaxation rate and leads to a Gaussian time-decay of spin polarization due to 
memory effects. In this case the relaxation becomes faster with increase of the magnetic field in 
contrast to the well known magnetic field suppression of spin relaxation. 



The effect of magnetic field on spin relaxation in low- 
dimensional structures in the presence of spin-orbit (SO) 
coupling is an interesting theoretical and experimental 
problem. The understanding of this effect which allows 
to a certain extent the engineering of the spin dynam- 
ics in these systems can have a crucial impact on the 
possible spintronics applications,^ The SO coupling of 
electrons in two-dimensional (2D) zincblende-based sys- 
tems with the structural asymmetry grown along the 
[001] direction is described by the Rashba Hamiltonian 
Hr = a{axky — cTykr^) (Ref. 3)i where a is the coupling 
constant, ai are the Pauli matrices, and k is the in-plane 
wavevector of the electron. The random spin precession 
necessary for the spin relaxation is introduced by elec- 
tron scattering by impurities, phonons, other electrons'^'^ 
and due to random dynamics in regular systems^ In the 
"dirty" limit, akr/Ti <^ 1 with t being the momen- 
tum relaxation time, the motional narrowing leads to 
the Dyakonov-Perel' mechanism^ causing the exponen- 
tial spin relaxation with the rate of the order of {ak/h)'^T. 
The spin relaxation rate decreases when a magnetic field 
is applied— due to the effect of the Lorenz force on the or- 
bital movement of the electron. From the analysis of the 
spin relaxation in a magnetic field valuable information 
both on spin and charge dynamics can be extracted^ In 
the clean limit, akr/Ti 3> 1 spin relaxation, on a time 
scale of the order of r, occurs on the top of the spin 
precessionSi. Various mechanisms of spin relaxation in 
solids are reviewed in Ref. 10]. 

In most quantum wells (QW) the Rashba-type SO cou- 
pling is achieved by asymmetric remote doping at the 
sides of the well; the same doping forms a smooth ran- 
dom potential scattering electrons. A crucially important 
step in being able to quickly manipulate the spins has 
been made by demonstrating that, by applying external 
bias across the quantum well, it is possible to change the 
magnitude of Q^iLi^iiiiiii^. Almost all previous stud- 
ies assumed the a parameter to be constant in space. 
However, the evidence is growing that the SO coupling 
both in zincblendc and Si^jGei-aj/Si-based QWsiSiii is a 
random function of coordinate. The randomness arises 
due to imperfections in the system, e.g. either due to 
shot noise accompanying the doping^^ or due to ran- 



dom variations of the bonds at Si/Ge interfacesii in the 
Sia;Gei_2,/Si systems considered as a hope for spintron- 
ics due to a very small SO coupling thereiiSiiS Therefore, 
the spin of a moving electron interacts with a randomly 
time-dependent effective SO field. This randomness leads 
to an irregular precession of spins and contributes to the 
spin relaxation rate. Here we concentrate on the systems 
where the Rashba term dominates in the SO-coupling, 
like InGaAs-based structures'^, and demonstrate new 
qualitative effects of the randomness on the spin relax- 
ation in a magnetic field such as a possibility of the non- 
exponential spin relaxation, and an increase of the spin 
relaxation rate in a strong magnetic field. The effect 
of the Dresselhaus coupling arising due to the unit cell 
asymmetry^'' will be briefly discussed in the context of 
our results also. The importance of the randomness of 
electric field of dopants causing the SO coupling on the 
donor electrons was first demonstrated by Mel'nikov and 
Rashba-i for bulk Si. Recently, Vagner et al.^^ showed 
that a random SO coupling for electrons in GaAs meso- 
scopic rings interacting with nuclear spins can lead to an 
Aharonov-Bohm-like effect. 

Since we are interested in the description of random 
spin and charge dynamics, we consider as an example 
a QW of the width w extended between —w/2 < z < 
w/2, and surrounded by two (5-doped layers at z = ±zo, 
zq 3> w/2, with a mean two-dimensional concentration 
of dopants (upper layer) and Nd (lower layer) with 
charge JeJ. The local SO coupling, being a linear response 
of the system to the symmetry-breaking perturbation, is 
presented as aR(p) = aso\e\Ez{p), where Ez{p) is the 
z— component of electric field of the dopant ions, p is the 
2D coordinate at the z = plane, and ago is a system- 
dependent phenomenological parameter. The z-and in- 
plane i = x,y components of the Coulomb field of the 
dopant ions are given bjtS^ 
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where Vj = {rji\,Zj), rj\\ = {rx,ry) is the 2D radius- 
vector of the jth dopant ion, zj — ±zq for the up- 
per and lower layer, respectively, e is the dielectric 
constant, Jo{qp) is the zero-order Bessel function, and 
Aq = exp(-gzo)/(9 + 9s), where qs = {2Tre'^ / e) dN / de f 
describes the Thomas-Fermi screening of the in-plane 
Coulomb field with TV = iV„ -|- Nd being the electron con- 
centration, and sp being the Fermi energy. As a result, 
electrons move in a smooth 2D random potential U{p) 
with -E||(p) = — V||J7(p), with the mean value {U) — 0. 

We consider the semiclassical movement of the electron 
in external magnetic field H: 
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where m is the electron effective mass. The momentum 
relaxation time for the electron moving in a smooth ran- 
dom potential is given by24 
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where Cjjuip) = {U{p)U{Q)) is the correlation func- 
tion of the random potential, and vp — hkp/m is the 
Fermi velocity corresponding to the Fermi wavevector 
kp — V2ttN. The relaxation time can be estimated 
as r TdE^p / lU'^) , the potential fluctuation ([/^) ^ 

N (e^/eg^^o)^, and is the time of passing through one 
domain of random force and random SO coupling, with 
typical size zq in this case. In InGaAs based structures 
with m ~ 0.02mo, where wlq is the bare electron mass, 
qs is close to 2 x 10^ cm^^. For numerical estimates we 
use the parameters from Refs. |13ii25<l with — 10^^ 
cm~2 and A'^d = 3 x 10" cmr"^. Then the ratio sj,/{U^) 
is of the order of 100 leading to r two orders of magni- 
tude longer than 0.03 ps. In weak fields lOcT <C 1, 
where uJc = \e\H/mc is the cyclotron frequency, the elec- 
tron path is random, while at lOcT > 1 it becomes close 
to the regular circular movement in magnetic field with 
a random path of the orbit center. At m = 0.02mo the 
crossover from the random to the regular movement oc- 
curs at H ^ 0.1 T. The semiclassical description of the 
orbital movement is possible in a non-quantizing field up 
to iJ - 1 T, where fiujc/ep ^ 0.1. 

The Rashba parameter is the sum of the mean value 
(a) = 27rasoe^(A^d ~ Nu)/e, and a random term with 
the zero mean contribution: a{p) = (a) + 6a{p) and the 
variations^ 
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S^{Nd-Nu)z, 
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The random and regular contributions become compara- 
ble at Nz"^ ^0.1; the condition satisfied in most of the 
experimentally investigated QWs. The regular term {a) 
can be reduced or completely removed by applying an 
electric bias across the quantum well, but even in this 



case the random term Sa{p) remains and causes spin re- 
laxation. Figure 1 presents a pattern of a^(p) obtained 
by a Monte-Carlo produced white-noise distribution of 
dopant ions corresponding to the experimental data of 
Ref. _13J. As can be seen in Fig. 1, the random varia- 
tions of QfR (p) are large and correlated on the distances 
of the order of zq. Therefore, the pattern of SO coupling 
is random nanosize domains rather than a constant. 




FIG. 1; (Color online) A realization of the random Rashba 
SO coupling (in arbitrary units) due to fluctuations of the 
dopant concentration, Nu = 10^^ cm~^, Nd = 3 x 10^^ cm~^, 
and zo = 12 nm. Curves (1) and (2) illustrate possible elec- 
tron trajectories without magnetic fleld and in the field in 
the ballistic regime, rc is the cyclotron radius and Id is the 
domain size (see Eq. 0). 

In general for a coordinate-dependent a{p), the 
Rashba Hamiltonian should be symmetrized to have the 
Hermitian form but we assume semiclassical movement 
of electrons and describe the spin dynamics with a mo- 
mentum and coordinate dependent effective Larmor fre- 
quency il[k{t),p{t)] =2a{p){ky,-ka:)/h. 

The semiclassical spin dynamics of electron neglecting 
the Larmor contribution of the external magnetic field^i 
to the spin precession is described by: 



ds 



sxn [k{t),p{t)] = 0. 



(5) 



Equation (jSJ means that the electron spin precesses in 
the effective field caused by the coordinate and wavevec- 
tor dependent spin-splitting that fiuctuates due to the 
variations of the wavevector and electron position. It is 
well known that the spin dynamics is governed by the 
fluctuations of the Larmor frequency CI provided that 
their timescale is shorter than the spin relaxation time. 
Let us introduce the autocorrelation function Cfio(i) = 
{ft[k{t),p{t)]n[{k{0),p{0)]) which is taken along the 
semiclassical path p{t) obtained by integrating Eq. 
Further we restrict ourself to the T-approximation2Si2S 
with r given by Eq. It allows us to present the so- 
lution of Eq. (jSJ for the ensemble mean value of spin 



(s,(i)) = (s,(0))exp 



- dt' dt"Cnn{t' -t") 



(6) 
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In high-mobility structures the smooth random poten- 
tial of dopant ions causes two main sources of random- 
ness leading to a spin relaxation: (i) random a (p) due 
to fluctuations of the electric field Ez and (ii) random 
direction of k. The former contribution causes the ran- 
domness in the amplitude of the SO field and, in turn, 
the precession rate, on the time scale r^, while the latter 
causes the random orientation of the precession axis on 
the much longer time scale r. These effects determine 
the correlator in Eq.©. 

To take into account different possible mechanisms 
of randomness of the SO coupling (fluctuations in the 
dopant concentrations, randomness of the bonds at the 
interfaces, variations of the directions of the growth axis, 
random strain, etc.) we consider a simple model with 



{da{p)Sa{0)) = {{da)^)e-p/^\ 



(7) 



where Id = r^vp is the characteristic domain size of the 
order of 10 nm. The details of the shape of the correlator 
in Eq.|(2} have no qualitative influence on our results. 
One can show that the correlator of the effective Larmor 
frequencies in Eq.® is given by 



n 



(8) 
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Here Sr is the displacement of the center of cyclotron 
orbit due to drift in the random potential, ^^'^^ and 
Ckk{t) = cos Wei e^*/'^ is the electron momentum corre- 
lator in the magnetic fleld. The first term in the square 
brackets gives the contribution to the spin dynamics due 
to regular SO coupling, the second one stands for the 
contribution from random fiuctuations of a. In Eq. (jSJ 
the correlations between the fluctuations of SO coupling 
and momentum are neglected since they have drastically 
different time scales and r, respectively. 

To begin we consider the case of a not very strong field 
^cTd ^ 1, where 5r ^ Id- In the collision dominated 
regime {ajkr/h <C 1, spin decays exponentially and the 
relaxation rate obtained with Eqs. © and Q is given 



by Tz 



-L ZZ 



where 
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-Td, 



{1 + UJ^T^)' 
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At _ff = both terms are the product of the typical spin- 
splitting squared and correlation time corresponding to 
the D'yakonov-Perel' relaxation mechanism.^ Here the 
electron passes many domains moving along the straight 
path (line (1) in Fig. the spin splitting in each domain 
is different and the time scale of the variations of the ef- 
fective Larmor field CI [k{t) , p{t)] is t^. In a magnetic 
field at w^T^ » 1, in the case of the regular SO coupling, 
the effective Larmor frequency S},[k{t + T/2)] — —i},[k{t)] 
is reversed each half of the cyclotron period on the cir- 
cular trajectory causing spin precession backwards, and 



the relaxation slows down. It is not the case when the 
SO constant fluctuates in space and (a) — 0: if the cy- 
clotron radius is large enough {ujcTd ^ 1), the random 
contributions to fl[k{t), p{t)] are not correlated, the re- 
versal of precession does not take place making the relax- 
ation rate iJ- independent. We emphasize that the source 
of randomness here is not the momentum scattering but 
the random fluctuations of the spin-orbit constant. A 
numerical estimate at {{Sa)^) = 1.5 x 10^^° eVcm and 

experimental conditions of Ref.fl^l gives l/vi'^z 60 ps 



To understand the qualitatively different situation 
when the magnetic field becomes so strong that after a 
cyclotron period electron returns to the same domain of 
spin-splitting, we consider the 'compensated' case with 
(a) = 0. For the particle in a smooth potential the con- 
dition for such a ballistic regime is {\p{t + T) — p{t)\) = 
5r <^ Id, T = 2i:/lOc, which can be achieved at cJcT 2> 
{t /TdY^^ This condition can be fulfilled in non- 
quantizing fields at T 3> 10 ps, corresponding to high 
mobilities /i ^ 10^ cm^/(Vs). In such a case on each cy- 
clotron period the electron passes exactly the same con- 
figuration 5a{p) (see circle (2) in Fig.^. It means that 
the effective correlation time of Larmor frequency would 
be much longer than and the spin relaxation becomes 
faster in a strong magnetic field. This result is in a sharp 
contrast to the well-known magnetic field suppression of 
the spin relaxation^ in Eq. (|5J| . During the first cyclotron 
revolution the spin dynamics is the same as in the ab- 
sence of a magnetic field. For subsequent revolutions 
the spin rotates by the same angle as during the first 
one because the electron passes the same configuration 
of a{p). This memory effect enhances the spin relax- 
ation rate on each cycle and makes the spin relaxation 
non-exponential, namely 



exp {-rW [K^T + (1 + 2K){t - KT)] } , 



(10) 

with K = \ t/T\ and [. . .J standing for the integer part. 
It is interesting to note that Cremers et alM found that 
the effect of non-uniform SO coupling on weak localiza- 
tion in quantum dots can be interpreted as a speeding 
up of spin relaxation. The cusps in {sz(t)) at i = nT are 
smoothed out if one takes into account the drift of the 
cyclotron orbit Sr in the random potential. At t 3> T 
we obtain {sz{t)) ~ exp (—Vi'^Jt'^/T). Figure 2 shows the 
calculated time dependence of the (s^) for different values 
of the external magnetic field. 

We note that in the high-mobility structures the col- 
lision dominated regime can be violated and spin relax- 
ation occurs on the top of spin precession making Eq. © 
invalid.^ However in strong magnetic field tOc ^ (ak/h) 
the applicability of Eq. 10 can be justified, besides, the 
regular contribution to the spin relaxation is strongly 
suppressed by the magnetic field. Therefore the spin 
relaxation rate drops to its residual value of the or- 
der of riz^ solely determined by the fluctuations of ran- 
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FIG. 2: Time dependence of (sz) in the compensated ballistic 
regime without magnetic field (solid), H — 0.1 T (dash) and 
H = 0.5 T (dash-dot). The arrows show the first cyclotron 
revolution. The parameters of calculation were as follows: 
m = 0.04mo, {{Say) = 1.5x10"^° eVcm, la = 12 nm, kp = 
1.6 X 10^ cm~^, and (sz(0)) = 1/2. This set of parameters 
can be realized in the structure investigated in Ref. H^ at 
iVa = A''d = 2 X 10" cm-^ 

dom SO-coupling. The further increase of magnetic 
field to LOcT (t/tii)'^^^ makes the spin relaxation non- 
exponential (see Ea. (|l()(l ) with the characteristic decay 

time Tg ^ \J T jV^z). ■ Finally, we mention that the ad- 
dition of the regular Dresselhaus SO coupling will not 



change our main results since it is suppressed by the mag- 
netic field in the same way as the Rashba contribution, 
and, therefore, only the effect of a random Rashba cou- 
pling remains in a sufficiently strong magnetic field. 

To conclude, we have shown that the inevitable ran- 
domness of SO coupling has an important effect on the 
spin dynamics in 2D structures, e.g. it determines the 
minimum spin relaxation rate achieved in a magnetic 
field. In the ballistic regime, achieved in sufhciently 
strong magnetic fields, the domain pattern of SO cou- 
pling leads to the non-exponential decay of spin polar- 
ization being a manifestation of a spin memory effect. 
The spin relaxation in a system with random SO-coupling 
becomes faster with the increase of the magnetic fieldi^ 
These effects limit possibilities of spin manipulation. The 
predicted magnetic field dependence of the spin relax- 
ation can be observable in 2D structures'^ with /i > 10^ 
cm^/(Vs), where the "ballistic" regime is achievable in 
non-quantizing fields. 
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